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The electromagnetic fields produced by an electric dipole immersed in an anisotropic medium are 
considered. Various approaches to the problem are outlined with special reference to a cold plasma. 
An attempt is made to show the close relationship between previously published work on this subject. 
It is shown that information on the radiation field in an anisotropic media may be obtained directly 
from the shape of the refractive index surface. 

1. Introduction 

There is now a vast amount of literature on the subject of electromagnetic waves in anisotropic 
media. However, until quite recently, only plane wave solutions have been considered in any 
detail. In most applications to the ionosphere, this approach has been justified since the trans- 
mitting and receiving antennas were ground based. In recent years, the situation is gradually 
changing since, increasingly often, the source is placed within the anisotropic medium. The best 
example is when a transmitter is located in a satellite which has an antenna attached to it in some 
manner. 

It is the purpose of this paper to review some of the recent theoretical investigations on the sub- 
ject of radiation from sources in anisotropic media. It is assumed throughout that nonlinear 
phenomena may be ignored. Thus, the medium is characterized by a tensor dielectric constant 
and a tensor magnetic permeability whose elements are independent of the field. 

In the first part of this paper, the Green's function for a generally anisotropic (i.e., gyrotropic) 
media is considered. The results are then specialized to a medium whose magnetic properties 
axe isotropic. This important special case corresponds to a cold plasma which may be described 
as a gyroelectric medium. Using this type of medium as an example, some alternative approaches 
to solving for the fields are described. The extensions required to treat stratified anisotropic 
media are also described briefly. Finally, the special problems associated with asymptotic evalu- 
ation of the integrals are considered. 

2. Green's Functions for Generally Anisotropic Medium 

The present development of the Green's function is closely related to that of Bunkin [1957], 
Kogelnik [1960], Meecham [1961], Kuehl [1962], and Ament [1963]. In particular, following 
Ament [1963], we assume a Cartesian coordinate system (xi 9 #2, #3) in which £}, Hj, and// represent 
the component of the electric field, magnetic field and source current density, parallel to the Xj 
axis. As usual, all field quantities vary according to exp {icot). 

The medium itself is considered to be homogeneous and generally anisotropic. Thus, the 
dielectric constant and magnetic permeability are denoted by tensors (e) and (/u,), independent of 
position. For this situation, Maxwell's equations for the jth component are written 

(VX~3)j + iG)fJL jn Hn=0, (1) 
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(VxH)j-io>e jn E n =Jj, 






(2) 



where fJLj n = (lJL) and ejn = (e) are three-by-three matrices or dyads with constant elements. 

Formally, the field E\[x) at the vector position x may be obtained via a Green's function 
Gij(x, x') from the current density Jj(x') at vector position x' . One should note that the compo- 
nents of x are Xi, X2, x% and those of x' are x[, x' 2 , x' 3 . The Green's function Gij is a tensor and is a 
function of both x and x' . It is defined by the relation 



Ei (?) = [ f f Gy (£ **) Jj (?) dx[dx' 2 dx'„ 



(3) 



where the integration is over any volume which encloses the source current density. A more 
compact description of (3) is written 



£«(*)= I Gf&lfijjbWx'. 



To obtain an expression for the Green's function, it is convenient to write Maxwell's equations 
in symbolic form. For example, if dt symbolizes the derivative djdxu it is evident that 



V XE=-i(jJL)(oH 



(4) 



is equivalent to 




£1" 


/Mn 


^12 


Mia\ 

M23 




"//l 


E 2 


= — id) I JIX21 


M22 




# 2 


.£3. 


\^31 


LL.Y2 


/W 




L#3 



or, simply 

d ■ E =— icojl • H. 

Similarly, the other of Maxwell's equations, 



(5) 



(6) 



(7) 



becomes 



d'H = ia)€'E + J. 



(8) 



The vector H may now be eliminated from (6) and (8) to obtain 

d • v - d • E — co 2 e - E = — icoj, 



(9) 



^here 



v = — is the matrix inverse to u. 
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To facilitate the solution of (9), it is assumed that both E and J have three-dimensional Fourier 
representations. These may be regarded as superpositions of plane waves propagating as exp 
| — ik\%\ — ikix-z — ik,ix.i] = exp [— ik • x] where k is the (vector) wave number. Thus, 



e n (Jc) = (2tt)- 3 I E n (x)e + * ' 7 d 3 x 



(10) 



is the transform of E n (x) when n= 1, 2, or 3. The inverse transform is given by 

En(x)=f e n (k)e-**(Pk, (11) 

which is a volume integral over. k space. The transform of J n (x\ which is denoted by j n (k), is 
defined in the same fashion. 

The differentiations in (9), indicated by the d tensors, may be carried out under the integral 
signs. Then, on equating coefficients of exp ( — ik • x), we obtain 

K • v • K • e + co 2 e • T= iojj, (12) 

where K represents the matrix 




Equation (12) may be solved formally to yield 

~e{k) = i(x>M T ](~k), (13) 



where 

M=~ —~ (14) 

K-v-K+aK 

From (11) we then find that 

~E(x) = ico J e~** • ~*M ^}(k)d 3 k ( 1 5) 

= m J e-M^M' J J(x r )e + ^^'{27T)- 3 dVd 3 k. (16) 

This result can be written in the required form 

~E(x > )=[ GU, ?)-7V)dV, (17) 

if 

C&T, "?) = -£*- f e~ Tk ■ <?- ^Md% (18) 

\2mY J 

and provided the order of the x' and the k integration may be interchanged. In matrix notation, 
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the Green's function is written 

Gy(?, 7') = ~ J c-* (? - ? >AM»Jb, ( 19 ) 

which involves the A: — dependent matrix A/*j. 

It may be easily shown that M-,j is an even function of k in that it has the same value when k 
(or hi, /c 2 , ks) is replaced by —A: (or by — A;i, — A; 2 , — fe). Thus, on noting that Afy is an explicit 
function of (x — x f ) and k, it is seen that 

djix- J?) = Cy(? - *) . (20) 

3. Green's Function for Gyroelectric Medium 

The development of the Green's function as given above is applicable to homogeneous but 
generally anisotropic media. The important special case is when the magnetic properties are 
isotropic. In this instance, 

Alo \ 

(A) = m* , (21) 

\0 /xo/ 

where fio is the usual magnetic permeability of the media. In most ionospheric-type plasma, it 
is permissible to regard /xo as the permeability of free space. However, this is not an essential 
restriction on the theory. 

A further simplification may be made if the d-c magnetic field is directed along the z axis of 
the Cartesian coordinate system. The tensor dielectric constant then has the form 

(22) 




where the (complex) constants e, e, and g depend on the properties of the medium. 

While the specific form of the Green's function for the gyroelectric case may be obtained as 
a special case, it is instructive to derive it directly following the method of Kogelnik [I960], 

The relevant wave equation, corresponding to (9), is now written in vector-dyadic notation, 

(vv-Al-^£)if=-;,xocj! (23) 

— > -> „ -> -* 

where A= v*V is the Laplacian operator, 1 is the unit matrix or dyad, and V V and e are also 

dyads. The elements of V V are the differential operators d 2 ldxid%k. 

Plane wave solutions of (23) have the form 

if(r) = J^e-^- ? , (24) 

Then, since 



dx, 
it follows that 



e -tk ■ x = _ ik . e -ik ■ 7 ; (25 ^ 



(v~V - Al-^j) e-' ? ' ? = - Q,Jk)e-^ ■ ? , 



(26) 
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where 



&(k) = tt-B+^- (27) 

C € 



is the "wave matrix." The homogeneous wave equation thus has the form 

nikWoe"* ■ * = 0. (28) 

The condition for nonvanishing fields is 

det tl(k) = 0, (29) 

where det denotes the determinant of the elements of the tensor. The planes of equal phase are 
specified by the wave normal n = (ni, ri2, ris) such that k = kn. The components of the direction 
cosines may be expressed in terms of angles a and (3 by 

rii ^sin a cos /3, 
7i2 = sin a sin /3, 
713 = cos ol. 

Using the specific form of the dielectric tensor defined above, it is not difficult to show that 

det h(t) = /i, G> 2 (€ sin 2 a + e cos 2 a)(k 2 - k*){k 2 - k\) = 0, (30) 



where 



where 



k i] = (^-g 2 ) sin 2 a + €€(1 +cos 2 a)±Y 

Af, J 2(c 2 /co 2 )(e sin 2 a + e cos 2 a)e ' (31) 



r = [(e 2 -^ 2 -€e) 2 sin 4 a + 4g 2 e 2 cos 2 a] * . (32) 

The roots, k x and A n , of (30), are the (complex) wave numbers for the "ordinary" and the "extraor- 
dinary waves. They depend on the angle a between the wave normal and the d-c magnetic field. 
Rather than using a Fourier representation to obtain the Green's function, it is instructive to 
utilize the properties of Dime's 8 functions. These are defined, in the usual manner, by 



I 



Ax')8(x-x')dx'=f(x) (33) 

J —00 

and 

8(x) = ^ J + %- ik ^dk 1 . (34) 

Then, if 

t (r)= [ GU , ?) '~7(x')<Px', (35) 
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it follows that G must satisfy 



V V -Al-^-H G(x,x') = -ifJLo(Dl8(x - %') , 



where 



8(x) = 8(x x )8(x 2 )8(x3). 



(36) 



(37) 



This step may be readily verified if both sides of the equation are multiplied by ,/ and then per- 
forming the volume integration. It should be noted that V V and A operate on the variables x 
and not x' . 

To complete the derivation, (26) is rewritten in the form 



V V- Al-^r ft-V- ifc •* = -le- < * •*, 



(38) 



where H _1 is the inverse of the "wave matrix." Both sides of this equation are now multiplied 
by exp {ik • x') and the integrations I d?k = I I dk\dk 2 dk 3 are performed. On identifying this 

equation with (36), it readily follows that 



<**,*) (27r)3 



\ n-'e-^^-x'Wk, 



(39) 



which establishes the desired result. As indicated, this is a special case of (19) for the generally 
anisotropic medium. Equation (39) was also given by Bunkin [1957] in slightly modified form. 
The explicit form for the inverse of (1 is obtained from Cramer's rule; thus 



n- 



A 



detft ' 



(40) 



where the matrix A is the "adjoint" of ft. For a gyroelectric medium of the kind described, it is 
found that 



A(k) = k 2 ~~k~k =k 2 ^ J+— 4 S\ 



where 6' is the adjoint of e , given explicitly by 

ee ige 

//€*= -ige ee 

e 2 



(41) 



(42) 



and J J e = 



' e(n\ + /x|) + e{n\ + n\) — ig(n\ + nl) + en x n 2 en x n 3 — ign 2 n 3 \ 

ig{n\ + rcf) + en^n 2 e{n\ + nf) + e{n\ + n|) en 2 n 3 + ign x n 3 

^e/ii7i 3 -h ign 2 n 3 en 2 n 3 — ign x n 3 e(l + n|) 



(43) 
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4. Modified Field Solution 

A three-dimensional transform technique for solving the dipole problem has also been given 
by Mittra and Deschamps [1963]. Their results have a number of unique features which should 
be mentioned. In particular, some ingenious manipulations are made in order to carry out the 
inversion of the transforms. Furthermore, their final expressions are not restricted to the far zone 
or to a lossless medium. 

The starting point, in Mittra's and Deschamps' analysis, is Maxwell's equations written for 
impressed electric current sources J and the tensor dielectric constant e as given in (22) while the 
permeability fjio is assumed to be that of free space. The eliminating H from Maxwell's equations 
leads readily to the usual inhomogeneous wave equation for E given by 

"V X V* X ~E - jjL a) 2 iJ? = - i/Jioojf. (44) 

A coordinate system is now introduced which diagonalizes the tensor e. The unit vectors in this 
system are denoted Vi, cX, and Us 9 and are defined by 

U l =2- l l' 2 (x + iy), & 2 = 2- 1 l 2 (x-iy), and U 3 =z, (45) 

in terms of the unit vectors i, y, and z in the conventional Cartesian system. Jt is understood 
that the d-c magnetic field remains oriented in the axial £ or Uz direction. An arbitrary vector F 
may be written in the [/-system as 

f = (j s f { + U2F2 + U3F3, (46) 

where 

Fi=F*- H% % F* = F- V u and F,=~F~U,. (47) 

In terms of the components F d , F y , and F z in the conventional Cartesian system, it is seen that 

F, = 2-V'\F. r - iF y \ F 2 = 2- l ' 2 (F x 4- iF v \ and F-, = F z . (48) 

The diagonal nature of the [/-system is immediately evident from the identity 

e = UiUi€i + UtUiei + ^ 3 J7 3 €3, (49) 

where 

ei + e + g, e 2 = €-g, and e 3 = c. 
The corresponding expression for the gradient operator in the [/-system is 

V = W, + M + W 3 , (50) 

where the derivative operators are defined by 

*- 2 -" , (s- i i)-*= 2 -"'(s + 'l>" d *-5- < 51 » 

After some algebraic manipulations, the inhomogeneous wave equation for E may be written in 
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convenient matrix form. Thus, in terms of the diagonalizing coordinate system, 



{-(duk + d® d\ dtds 

d\ -(rfiA + dJ) d 2 d s 

d 3 d 2 did* —2d\dii 



-/x c*r 



€1 





°\ 


Ex 




Jl 





€2 





E 2 




h 








*/ 


E 3 




L 



(52) 



These are matrix equations for the fields,' first given in this form by Mittra and Deschamps [1963]. 
To solve these, it is necessary to employ Fourier representations which transform the differential 
operators d\,d 2 , and d$ to algebraic quantities. 

The three-dimensional Fourier transform of £i(i = l, 2, 3), denoted #ii= 1, 2, 3), is defined by 

#i = (2tt)- 3 "I J + °° Ei (*, y, z) exp [+ i{k x x + k y y + k z z ] dxdydz . (53) 

J— oo J— oo J— 00 

The differential operators now transform in the following manner 

rfi -> 2- 1 / 2 i(A: J . - iky), d 2 -» 2- l >H{k x + ik y \ d 3 -> ife. 

For convenience of field calculations, it is desirable to introduce the usual polar form defined by 
kx — k sin a cos /3, k y = k sin a sin /3, and k z — k cos a, whereupon 

rfi — » i2~ ll2 k (sin a)e"^, ck ~ > i2~ y/2 k (sin a)e^, fife --* i/c cos a. 
The transform of the matrix equation (52) now has the form 




1 



|f A: 2 1 ^¥^ 



s 2 e- 2 ^ 



f + c* 



v2- 1 ' 2 sce + # 2- 1 ' 2 sce-^ 




\6\~ 




'■Si 


► ^ 2 


= — ifJLoCt) 


<f-l 


Us. 




A. 



, (54) 



where s = sin a and c = cos a. 

The determinant A of the coefficients inside the curly brackets of (54) is seen to be of the form 



A = /x, co 2 (e sin 2 a -he cos 2 a)(k 2 — k 2 )(k 2 — A; 2 ), 



(55) 



which is identical to (30) as it must be. The roots of the equation A = are k x and k u and they are 
given explicitly by (31). As indicated previously, kjko and k n /ko are the indices of refraction for 
plane ordinary and extraordinary waves, respectively. 

Since the matrix equation (54) involves only algebraic quantities the solution may be obtained 
directly. Thus, in matrix notation, 



[f]=-i^co[A][jfl 
where the elements of the inverted matrix are given by 

A rnn ~ & mnl A • 

The coefficients a mn are given by 



an = (1/2)AV - k 2 fjLo<*) 2 \ e 2 s 2 + e 3 ( - + c 2 ) } + /x co 4 e 3 e2, 
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(56) 



(57) 



fli2 = e-W{\l2WsW ~ /x co 2 € 3 ), 
ai 3 = e-^2- l l 2 k 2 sc{k 2 - p oj 2 e 2 ), 
a 21 = e 2 ^(l/2)Fs 2 (£ 2 - /*oco 2 e 3 ), 



a 2 2 = (ll2)k 4 s 2 -k 2 fjL co 2 



+ e 3 (tt + c 2 



-hjLt 2 w 4 eie 3 , 



a 23 = e^2-^ 2 k 2 sc{k 2 - /lkk» 2 €i), 

a 3 l = e i n-^ 2 k 2 SC{k 2 - /JLo(D 2 €2), 

a 32 = 2-^ 2 e-^k 2 sc{k 2 - )Ltoa> 2 ei), 

a 33 = & 4 c 2 — /ULo(x) 2 k 2 (e 2 + d) (J- + c 2 J + /x|a> 4 €i€ 2 . 

The explicit evaluation of the three-dimensional inverse transforms remains to be done. The 
usual approach is to carry out a stationary phase type calculation which is stated to be valid in the 
far zone. This aspect of the subject is discussed below, but first, it is appropriate to describe the 
approach by Mittra and Deschamps [1963], which is not restricted to the far field. 

To be explicit, an electric dipole source is located at the origin and, to achieve simplicity, the 
dipole is oriented in the z-direction which is parallel to the d-c magnetic field. Thus 

~T= Po z8(x)8(y)8(z\ (58) 

where z is a unit vector in the z (or (/ 3 ) direction. Thus, from (55), (56), and (57), 

/> ■ (WA;, a, ft) r 

6 n = — ilJL {) a)po — 7T} : — lor /i=l, 2, 3, (59) 

and the inverse transform may be written in the form 

E n =Po | \ | Sne- ikR[sin e sin a cos ( *-* )+ cos e cos °'X k 2 sin adfidadk, (60) 

Jo Jo Jo 

where (/?, 0, cf>) are the spherical coordinates of the observer. 

In the special case n= 1, the /3 integration is effected immediately since a 33 is independent of /3. 
Thus, it follows that 

where 

/(/?, 6) = ^ 2 | | ^p- e"*** cos * cos «Jo(A:/? sin sin a) sin adadk. (62) 

Jo Jo A 

To evaluate this integral, Mittra and Deschamps [1963] recognize that / may be written as a sum of 
a singular part I s and a finite part //. That is, I s is of order (l/R) for R tending to zero, whereas 
// is finite for all R and 6. Following some ingenious manipulations, they find that 

T72r * M 1/2 l d d l 

V2/ < = ~2l7j ~ P Tp P TpW> ™ 

127 



where R f : 



e 
— P* 



72 

This result has the correct singularity and, in addition, it demon- 



strates that I s is independent of g, the off-diagonal factor in the expression for the tensor dielectric 
constant (e.g., (22)). 

The evaluation of V 2 // appears to be more complicated. It turns out that while // is finite at the 
origin, \Hf does have a singularity of the type 0{l/R). Mittra and Deschamps [1963] have given an 
explicit expression for V 2 I f but it involves double integrals which, fortunately, have finite limits. 
Apparently, a great deal of numerical work is required in order to obtain useful numerical results 
from this approach. It is rather surprising that they do not carry out the ct integration explicitly 
by using contour integration methods. This requires that one evaluate the residues of the poles 
of the integrand which occur at the zeroes of A for a fixed value of k. This operation will not 
be carried through since it leads directly to the (single) integral representation, given by Barsukov 
[1959] and Arbel and Felsen [1963], which may be obtained in a more direct fashion as outlined 
below. The latter approach is particularly appropriate when dealing with horizontally stratified 
media. 

5. Sources in Stratified Anisotropic Media 

Before tackling the generally stratified medium, the appropriate expressions for the (primary) 
field of a point source in a homogeneous medium are derived. With respect to the Cartesian co- 
ordinate system, an electric dipole is located at (0, 0, z ). The d-c magnetic field is oriented along 
the z axis and thus the dielectric dyadic e has the tensor form. 



(64) 



where e, e, and g are specified constants for each homogeneous region. For convenience, Max- 
well's equations 

~Vx~H = ia)€~E+~7 (65) 



V X E =-ifjL o)H (66) 

are supplemented by the divergence relations 

V(e- £) = p = — ~V •? (67) 

(O 

and 

V •#=(), (68) 

where p is the charge density. 

The preceding equations may be written out explicitly in terms of cylindrical coordinate sys- 
tems. Then, it is a straightforward matter to eliminate the E r , E^, H r , and H^ components. It 
is found that 
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and 



vmer 



and 



e d 2 E z KkiE z fJWg BH Z . . 1 dp 

V , H-+ ^ + 5. (e ._^ ff . + i^— ^x?]. + ^ , (70) 

dz 1 e<Eo e dz e 

z dp p d(f> 

V 2 = IA p ^ + lJ_ ( 72) 

p dp dp p 2 d 2 (/> 

is the transverse Laplacian operator. In this section k = o)/c = 27r/wavelength and e is the di- 
electric constant of free space. 

The above coupled equations for E z and H z will hold for any (piecewise) homogeneous region. 
A subscript m is added to the field quantities when reference is to be made to the rath homogeneous 
region. The integer m is allowed to range from 1 to M. The source dipole is taken to be in the 
pth layer where p is an integer equal to or between 1 and M. 

To facilitate matching of boundary conditions, it is convenient to employ the following spectral 
representations of Dirac delta functions: 

1 f + ac 
8(z-2o)=— e iu(z ~ z ^du (73) 

277J-3C 

^ £1 =\ Mkp)kdk. (74) 

P Jo 

Therefore, if the dipole of moment p , located at z = z , is oriented in the z direction the current 
density has the following form: 

A = /:a>poSU-2„)8(p)/(27rp) = -^ f* [^ J (Xp)e"«*-^kudkdu. (75) 

(Z77)~ Jo J- x 

The corresponding formula for the charge density is 

P = ~-^b " I + " Jo{kp)e***-*d\ud\du. (76 ) 

\ZiTT) J o J - x 

In the sth region, which contains the source, the axial field components are written in the form 

E,=Ep+E° (77) 

and 

H z =Hf+H>, (78) 

where the primary and secondary fields are indicated by a superscript p and s, respectively. It 
is understood that, if the 5th region were unbounded, the secondary fields would not be present. 
Taking a clue from the appearance of (75) and (76), the axial components of the primary field 
are written in the form 



and 



£? = 7^i f " f + ~ A{\,u)J () {\p)e h ^-^d\du (79) 

ff?=7^T s f °° f ^ B(K u)j»(kp)e"^-^d\du. (80) 

(Z7T) Z J o J -oc 
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Substituting (75), (76), (79), and (80) into (69) and (70), explicit expressions may be found for the 
functions A and B. The integrations, with respect to u, may be easily carried out and the final 
result is 



and 



^her 



6 Jo 5 2( W 2 u \) 



Hp=-*£ Dn (z - Zo) ( c °\e-»i\*-*o\-e-<'2\"-*ol*\X Uk ( ,)X3dK 



(81) 



I [ e -«i|«-«ol — e~ M 2l*-*ol]x 

Jo 



P07 ^ [ e -"»l*-*ol- e -tt2l*-*ol]x , 9 , (82) 



I z — z 1 J o u\ — u\ 



S 2 = ^2 (e/eo) _ X 2 ? 5 2 = (e/ ^) [*2(€/€ ) ~ X 2 ] , (83) 

where Wi and a* are roots of the equation 

JAo-2 
u 4 + n 2 (s 2 + 5?) + S 2 S 2 - 5 ? J2£_ = o (84) 

€6 

chosen so that Re u\ > and Re uo > 0. 

Expressions for the other components of the primary field may be readily found from the 
axial components given here. 

Expressions for the secondary field, in the 5th region, are now sought in the form 



E s z = ^ z Po f 3C {wi(a?-h5?)[^i(\)e- M ^ + i/i 1 (X)e +M ^] 

6 Jo 



+ wAiii + s 2 ) |> 2 (\)e - "** + i// 2 (X)e + wj} ^f^ (85) 



and 



H z = -^fp f "{[^(XJe-^- WX)e +MlZ ] - [> 2 (X)e-"2*- W x)€+V]}^^^, (86) 

where the functions <£>i(V), ^(A-), i//i(X), and ijMM are undetermined. 

The axial field components, E n z and //", in the other homogeneous slabs would have the same 
form as the above expressions if the appropriate material constants of the medium are employed 
(e.g., for the nth slab e, e, and g are replaced by e n , € n , and g n , respectively). In each slab, there 
are four unknown functions (i.e., <£>?, (p%, i/jjf, and *//??), and these may be found from the matching 
of the four tangential field components at each interface. In addition, it is required that only 
outgoing waves are permitted at z = ±°°. Thus, in principle, the complete formal solution is 
available. 

Using the method outlined here, Barsukov [1959] has given an explicit expression for the 
field of a vertical electrical dipole at z = z , which is immersed in a semi-infinite anisotropic medium 
(z > 0). In the region below (i.e., z < 0), the medium is assumed to be free space with electrical 
constants e and /x . 

For z > 0, the axial fields have the form 

E t =E* + & g , (87) 

H z = Hp + H s z . (88) 

Here, E p z and H p z are the primary fields given by (81) and (82), while E s z and H s z are the secondary 
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fields which have the form of (85) and (86) with the condition that i//i(X) and i// 2 (X) are both zero. 
For the isotropic region z < 0, it is convenient to represent the field in terms of electric and 
magnetic Hertz vectors, Yi e and II m . Thus 



E = Argil* + V(VFU) - ifJLoioV x n m , 
H = ie Q o) VxII e + k$\ m + ie o) V X Yl e , 



(89) 
(90) 



^here 



rie=(o,o, lie) and n w = (o,o, n m 

The following representations are employed: 

n e =^ f" ' F 1 (k)J (kp)e^kdk 
«o Jo 



and 



^o Jo 



F 2 (\)J (\p)e u " z \dk, 



(91) 



(92) 



where i^o = (X 2 — A:©) 1 ^ 2 with Re u > 0. The scalar functions, Yl e and n m , satisfy the wave equation 
( V 2 -h k$)He,m — 0; they behave as outgoing waves at z = — °°, and they are independent of the azi- 
muthal coordinate 0. The functions <pi(A), <p 2 (M, Fi(X), and F 2 (X) can be determined explicitly 
from the four boundary conditions (requiring the continuity of the tangential fields) at z = 0. 
Barsukov [1959] gives 
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(pl(Uu U'l) 
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e8(ai, ui) 



U\(u\ + s'f)(^o + Ui)e~ u \ z o — w>(ui -h sj)(uo + ai)e~ w 2*o 
u\{eu\Ui — €os|)e~ w i*o — uzieuiUz — 6osl)e~ u 2 z o 



^here 



5(ai, ^2) — (u> — Hi) [e^sliuxWi — s\) — euliuiWi — s'i)-\- uq(ui -f uo^eos'i — eu]Uo)] 
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(94) 

(95) 
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This constitutes the complete formal solution of the problem. Explicit expressions for the 
field components are obtained by differentiation according to (89), (90), (91), and (92). 

The general results are rather complicated; however, some simplicity is achieved by consid- 
ering the radiation fields in the lower half space. Physically, this corresponds to the situation 
where the observer is at a distance, below the interface, which is large compared with the wave- 
length. The integrals to contend with are of the form 



/: 



A(k)J (\p)e- u Mkd\ = I(p, z), 



(98) 



where A(K) is a slowly varying function of X. A saddle point evaluation of the integral yields 



733-195 0-64— 3 



/(p, z) ~ ik cos 6A(k sin 0)e -**«*/#, 
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(99) 



where 

R = (p 2 + z 2 y' 2 and tan 6 = -{plz). 

As indicated, the saddle point or principal value occurs at k = k sin where 6 is the angle 
subtended by the direction to the observer and the negative z-axis. This asymptotic result is valid 
provided \koz\ > > 1 and |z| > > Zq. The first inequality appears to be very important since if 
kz is small, the contributions from the poles of A(\) lead to surface waves which have a significant 
amplitude compared with the space wave. The pole wave contributions were not considered by 
Barsukov [1959]. 

6. Radiation Fields in Anisotropic Media 

An important aspect of this subject is the nature of the radiation or far-zone fields when both 
source and observer are within the anisotropic medium. This particular problem has been con- 
sidered in some detail by Arbel and Felsen [1963]. In their work they make extensive use of the 
geometrical properties of the refractive index surfaces. For purposes of this discussion, the me- 
dium will be considered homogeneous and unbounded. Furthermore, the source is an electric 
dipole of moment p located at the origin of a cylindrical coordinate system (p, </), z) and oriented 
in the z direction. The applied d-c magnetic field is also directed in the z direction. The medium 
is nonmagnetic with permeability /x , while the electric properties are described in terms of a ten- 
sor dielectric constant e of the form given by (22). Losses in the medium are assumed to be zero. 
The Fourier-Bessel representations of the fields given by Arbel and Felsen [1963] are closely re- 
lated to (81) and (82) first given by Barsukov [1959]. The essential modifications to (81) and (82) 
are to replace the variable A by k&i where S is now a dimensionless variable and k = ojlc. By 
using identities Jo(x) = i [H$\x) H- H$\x)] and H^\x) = — H^\ — x), the integration contour may be 
changed from along the positive real axis to a straight line running from S = — °° to -f-°°. Because 
of a branch point singularity at S = 0, the integration path is indented around this point by a small 
semi-circle in the lower half plane of the complex S plane. Additional substitutions in (81) and (82) 
are to replace the axial propagation constants u, ui, and ui by ik C, HcqCi, and ikoC-z, respectively. 
The subscripts 1 and 2 may be identified with the ordinary and the extraordinary waves in the clas- 
sical magneto-ionic theory. 

The resulting vector electric fields may be written in the form 

~E=~E^+~E 2 

where the field components are given by 

(£ P )i, 2 = iA A \ +X S 2 ^^ B\\koSp) exp [- ikoCt, 2 \z\] dS, (100) 

\Z\ J— oo ^ 

(E <t> ) l , 2 = 2A~^ j + ~ — WfikoSp) [-ikoC h2 \z\] dS, (101) 

(£*)i, z=A J* ^^ y^ Hty&p) exp [- ikoCi, 2 \z\] dS, (102) 

A = ix co 2 poI(16tt), (103) 



and 



/here 



and 
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The functions C\ and Ci are proportional to ii\ and ui as noted. The latter may be obtained from 
the roots of the quadratic equation in ir given by (84), which is equivalent to the Appleton- 
Hartree equation. Explicitly, 
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Arbel and Felsen [1963] have devoted a great deal of attention to the proper definition of the 
multivalued function C and the disposition of the integration path in relation to the singularities 
of the integrands in (100), (101), and (102). In the first place, it is to be observed that C\ and C-i 
are associated respectively with the + and — signs in (105), and that {U 2 — W) 112 is defined to be posi- 
tive when real. The definition of Ci given above implies that C\ has branch-point singularities at 
U 1 = W (or at A = 0), and at those values of S for which Ci=0. Similarly, Co has branch-point 
singularities at A = 0, and at those values of S for which C 2 = 0. However, a series expansion of 
the total integrand about A = contains only even powers of A, which means that A = is a regular- 
point. Thus, no special care need be taken in the definition of Ci, 2 at A = if the "ordinary" and 
"extraordinary" are combined and treated together. Also, it should be mentioned that the factor 
(I — S 2 )~ l in the integrand of (102) does not give rise to a pole contribution since it is canceled by 
other factors. 

If C\ or C 1 are real, the quantity exp ( — ikoCi, >>z) represents a propagating plane wave in the 
positive 2-direction away from the source. On physical grounds, it may be supposed that each 
component plane wave, in the total spectrum, has a positive component in the z-direction. This 
radiation condition permits the unique determination of Ci,2 over the range of S for which the 
propagation constants are real. As Arbel and Felsen [1963] show, the contributions to the power 
flow, when the propagations are complex, are zero. The resulting "energy radiation condition" 
is that, when Ci, 2 is real, 



Mo 
60 
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(A 2 ±S 2 A) \€-l 



>0, 



(108) 



where Y\ and Y> may be interpreted as admittances for ordinary and extraordinary waves traveling 
in the positive z-direction. 

It is evident from the above inequality that the sign of C\ t 2 is identical with that of A 2 — S 2 A. 
This condition for real &, 2, when supplemented with the condition Im &, 2 < for C\, 2 complex, 
may be used to determine the analytic continuation of the function Ci,2<S) around branch-point 
singularities on the real axis of the S plane. 

As indicated above, the propagating modes correspond to the situation where 



Cl 2 = U±(U 2 -W) 1 ' 2 



(109) 



is positive real. Thus, it is evident that ordinary and extraordinary modes propagate when U > 0. 
On the other hand, if (U 2 — JV)> |£/|, i.e., W<0, the ordinary mode propagates for U^ 0, while 
the extraordinary mode does not propagate. 

There is actually a very simple graphical procedure which utilizes the dispersion curves of 
C{S) versus S. For example, the plane wave, characterized by the variation 



exp [ — ik (Sp 4- Cz)\ 
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carries power in the direction perpendicular to the dispersion curve at the point S, C. The radia- 
tion condition corresponding to choosing those values of C\ and Co for which the outward normal 
to the dispersion curve, has a positive component in the z-direction. 

The integrals in (100), (101), and (102) maybe evaluated asymptotically by the method of saddle 
points. This has been demonstrated by Arbel and Felsen [1963] in an elegant manner. 

It is convenient to introduce spherical coordinates (r, 0) which are defined by p = r sin and 
z = r cos 0. Then, under the assumption that k r sin 6 > > 1, the asymptotic approximation 

H}; 2 \koSp) ~ ( , c 2 , Y " exp [ - ikoSr sin + inirl'l + irr/i] (1 10) 

\7rkobr sin 0/ 

is introduced. Thus, for most practical purposes, the integrals to contend with are of the form 

B= \ L(S)exp[-ik rM(S)]dS, (111) 

where 

M(S) = C(S)\cos $\+S sin 6, (112) 

and k r is to be regarded as a large parameter. 

The principal contributions to the integrand occur at the vicinity of the saddle points S* deter- 
mined by 

M'(S*) = 0, (113) 



C'(S'') = -|tan0|. (114) 

Propagating waves correspond to the saddle points where both CiS 1 ) andS' are real. If these quan- 
tities are complex, the associated fields are highly attenuated and need not concern us further. 

The technique to evaluate the integral asymptotically consists of deforming the integration 
path to the steepest descent paths through each saddle point. The ends of these paths are con- 
nected up in the valley regions of the S plane where — Im C(S) is a large positive quantity. Thus, 
because of the heavy exponential damping, these portions of the path contribute a negligible 
amount and may be ignored for the present analysis. In the vicinity of a real saddle point S,, the 
factor exp [ — ik rC(S)] decays most rapidly along the steepest descent path which is inclined at 
angles of ±45° with the real axis. 

The asymptotic representation for the integral is then found to be of the form 

B = ]£ B(i) ' (115) 

where the contribution from an individual saddle point is 

2tt 



S»> ~ US'] 
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e +iir/4 e -i* rW*) 9 (H5) 



where the a r gument — J7r/4 is used when M"(S') > 0, and + i7r/4 is used when M'XS 1 ) < 0. It is 
assumed that the integrals are of the form such that the saddle points have alternating signs along 
the real axis of S. 

To express the final results for the radiation fields, the unit vectors, p, (/>, and i are introduced 
along with the polarization vectors 

]? = pS i +zO and^pff + f 1 "^ 2 • (117) 
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Then, following Arbel and Felsen [1963], it is found that the asymptotic versions of (100), (101), 
and (102) may be written 

„ , p-ikorNi(O) — » * 

pE p + zE z ~ 2iA X /3W /2 G f j—iq X 4>), (118) 

and 

E^2A^ p<JWF* e —j—, (119) 

where Ni(6) is the "ray refractive index" defined by 

MiS^NHO). (120) 

The radius of curvature R l of the dispersion curve (i.e., C(S) versus S for real S) at the saddle point 
S l is obtained from 

Rl = {l 7w)f }3/2 = {|r(S/) ' CQs3 e} ' K < 121 ) 

In (118) and (119) j3' = 1 if M"(S') < and j8* = i if M"(#) > 0. Finally, the amplitude factors are 
obtained from 



and 
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(124) 



It is understood that (118) and (119) both hold individually for ordinary and extraordinary waves. 
The summations extend over all contributing saddle points in each case. The exponential factors 
involve the ray refractive index N l {6) which is related to the polarization vector pt=pS i +zC' by the 
scalar identity N i {6)=p • r where r is a unit vector parallel to the ray direction which points from 
the source to the observer. 

For the results given above, it is specifically assumed that e/€ > 0. If e/e < 0, then the rele- 
vant field expressions are obtained by making the substitutions 

(eyl 2 ^-i(-e) 112 , C^iC, and S^iS. 

The characteristics of the radiation field in an anisotropic medium are critically dependent on 
the shape of the dispersion curve (or refractive index surface). It is evident that the number of 
contributing rays and their amplitude may be determined by a careful examination of the dispersion 
curve. The presence of open branches on the dispersion curve indicates presence of shadows 
where no rays may propagate. However, the field in the shadow region may be finite as a result of 
diffraction [Brekhovskikh, I960]. This aspect of the subject is outside the scope of the present 
survey. Also, the presence of points of inflection in the dispersion curve (i.e., C'XS 1 ) = 0) will lead 
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to focusing effects which may be described by a modified saddle point procedure involving Airy 
functions [Wait, 1962]. This topic is also outside the scope of this discussion. However, the 
interested reader may find an elegant treatment of this particular problem in the paper by Arbel 
and Felsen [1963]. 

7. Final Remarks 

The present expository paper has been a self-contained treatment of the theory of radiation 
from sources in cold magnetoplasma media. 

In this review, no mention has been made of antenna impedance for anisotropic media. While 
this would seem to be a natural topic for this discussion, it has been deferred until a later time. 
One reason for this is the need to consider thermal effects and sheath phenomena, which are 
important near the source. Consideration resulting from the finite temperature leads to contribu- 
tions from the acoustic-type waves not usually considered part of magneto-ionic theory. 



I am grateful to Dr. K. Furutsu for making a number of corrections in the manuscript. I also 
thank Mrs. Eileen Brackett for help in preparing this material. 
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